BMC Ecology and Evolution

Conspicuous coloration in animals is generally thought to evolve and be maintained through inter- or intraspecific interactions such as mate choice, but this might not always be the case. The sight-line hypothesis proposes that conspicuous light-dark contrast in front of the eyes (hereafter, eyeline) evolves and is maintained due to viability selection, enhancing an individual visual acuity and thus evolutionarily associated with a particular foraging behavior that requires accurate aiming. However, empirical evidence that supports the sight-line hypothesis is virtually absent, with no studies demonstrating the key prediction that the direction of eyelines matters. Here, I tested the sight-line hypothesis using macroevolutionary analyses in terns and allies, which are a suitable study system, because they have variation in facial color patterns, including presence/absence and, if any, various angles of eyelines. They also have a large variation in foraging behavior, including picking, plunge diving, and skimming. As predicted by the sight-line hypothesis, tern lineages that require accurate aiming at foraging (e.g., plunge diving) are more likely to have eyelines. In addition, the evolutionary transition to the state with eyelines and these foraging behaviors was more likely to occur than the reverse transition. Furthermore, as expected by the fact that the direction of travel is upwardly deviated from the direction of the bills during skimming, the eyeline angle from bills was evolutionarily positively associated with the occurrence of skimming behavior. To my knowledge, the current study is the first to demonstrate that the direction of the eyeline matters, thereby strongly supporting the sight-line hypothesis.

High and low tides occur twice a day (every [~]12.4 hours), with the largest tidal ranges during spring tides around new and full moons (every [~]14.765 days). While these lunar cycles are known to influence many animal phenotypes, particularly the reproduction of coastal animals, the genetic basis of lunar-related rhythms remains unclear. Since phenotypic variation is a valuable resource for elucidating such mechanisms, we examined geographic variation in the lunar-regulated mass spawning of the grass puffer (Takifugu alboplumbeus) along the Japanese coast. We found that western populations spawn during the first half of the spring tides, whereas eastern populations spawn during the second half. Furthermore, although spawning typically occurs a few hours before high tide, this timing is restricted to a specific time window that is earlier in the western populations than in the eastern ones. Behavioral analysis of larvae also revealed a shorter free-running circadian period ({tau}) in the western population than in the eastern ones. As differences in {tau} affect individual variation in the timing of physiological functions and behaviors, we hypothesized that differences in {tau} could account for the different time windows and consequently the observed difference in spawning days. Population genomics analysis identified proline-rich transmembrane protein 1-like (prrt1l) as a candidate gene. Expression of prrt1l was observed in the circadian pacemaker suprachiasmatic nucleus, and triple CRISPR F0 knockout of prrt1l shortened the free-running period in larvae. These findings suggest a potential mechanism underlying the geographic variation in lunar-synchronized spawning behavior. HighlightsO_LIThe geographic variation exists in the lunar-regulated spawning of the grass puffer, with differences in spawning dates and times between western and eastern Japan. C_LIO_LIThe free-running period of western populations is shorter than that of eastern populations, which is consistent with their earlier spawning timing. C_LIO_LIPopulation genomics analysis identified prrt1l as a candidate gene harboring population-specific missense mutations, the knockout of which shortens the free-running period. C_LI

Honeypot ants represent an example of convergent evolution, where a group of workers specialized in storing liquid food in their crops (i.e., stomach) has independently evolved multiple times across different ant genera. While seasonal resource scarcity and arid conditions are thought to drive the evolution of repletism, the role of environmental variables in this process has not been tested. With this is mind, species ensemble models were computed to assess suitability and richness areas, and the importance of predictors. Predictor importance was compared between genera and groups occupying a similar geographical area. Niche overlap and similarity between honeypot ant species were also evaluated to determine whether they occupy similar environmental spaces. Similarity was mainly found within genera, and Leptomyrmex and Myrmecocystus showed striking niche differences. Overall, Leptomyrmex distribution was mainly influenced by atmospheric bioclimatic variables like precipitation and temperature, while Myrmecocystus had soil bioclimatic variables as the most important predictors for their current distribution. Our results indicate that honeypot ants species currently do not occupy the same environmental space, and are not experiencing the same contemporary environmental stressors. While our results suggest that contemporary environmental factors cannot explain the convergence of honeypot ants, future research will examine past climatic conditions along with investigations into the ant genomes to understand more about the causes and consequences of the convergence.

Sarita Lake, British Columbia houses a distinctive population of threespine stickleback (Gastrosteus aculeatus L.) with a phenotype characterized by unusually large individuals relative to nearby conspecifics. We tested the hypothesis that members of this population are not isometrically larger but rather exhibit variation in allometric trajectories that reflect changes in developmental timing impacting the developmental-genetic architecture of the phenotype. We used 3D geometric morphometrics to characterize the size and shape of skulls, pectoral girdles and pelvic girdles from a sample of individuals from nearby freshwater and marine populations and compare them to a sample from Sarita Lake. We showed that individuals from the Sarita Lake population are larger in each body region compared to most other populations examined. Further, these individuals have dorsally expanded skulls and relatively robust pelvic armour. We also showed that the relationship between size and shape is differently structured among body regions and is heavily influenced by non-uniform sexually-mediated variation across populations sampled. Our results reflect complex underlying developmental trajectories, and we suggest that the large phenotype observed may be driven by fecundity selection on female size in combination with a limnetic trophic niche and relatively increased predation pressure in Sarita Lake.

Text AbstractIn this study, we characterized the genetic structure and reconstructed the demographic history of cactus wrens (Campylorhynchus brunneicapillus), an endemic species of desert regions of North America, that shows a clear phenotypic and genotypic variation. We evaluated the effects of historical climate change on the structure and population dynamics of desert species using genomic data through genotyping by sequencing (GBS) and applied a population structure analysis (FST and ADMIXTURE), revealing two genetically differentiated groups: one continental and another peninsular in Baja California. Subsequently, we implemented the MSMC2 coalescent model on data divided into autosomal regions and the Z sex chromosome to estimate changes in effective population size (Ne) through evolutionary time. Additionally, we developed ecological niche models (ENMs) projected to the Last Glacial Maximum (LGM), Last Interglacial (LIG), Present times, and Future (2060 - 2080). Results indicate that both populations maintained moderated Nes before the LGM, experienced severe bottlenecks (Ne [~] 102-103), followed by a sustained expansion. However, recovery was limited to the Z chromosome of the peninsular population. These findings reveal how glaciations and interglacials shaped the evolutionary history of desert species and provide genomic evidence of the splitting of C. affinis from C. brunneicapillus. Article summaryThis research examines how climate changes shaped genetic diversity of cactus wrens across North American warm deserts. Using coalescent methods, researchers tracked effective population size changes over 100,000 years, using ecological niche modeling they predicted habitat suitability across climate periods. Results showed that continental and peninsular populations experienced bottlenecks during the Last Glacial Maximum, followed by demographic recovery on warm periods. However, the sex chromosome (Z) revealed male-biased demographic patterns in peninsular populations. Future projections indicated habitat suitability reductions for peninsular populations, highlighting conservation concerns. These findings demonstrate that past climate shaped genetic diversity of cactus wrens.

At the fundamental conceptual level, two alternatives have traditionally been considered for how mutations arise and how evolution happens: 1) random mutation and natural selection, and 2) Lamarckism. Recently, the theory of Interaction-based Evolution (IBE) has been proposed, according to which mutations are neither random nor Lamarckian, but are influenced by information accumulating internally in the genome over generations. Based on the estimation-of-distribution algorithms framework, we present a simulation model that demonstrates nonrandom, non-Lamarckian mutation concretely while capturing indirectly several aspects of IBE: selection, recombination, and nonrandom, non-Lamarckian mutation interact in a complementary fashion; evolution is driven by the interaction of parsimony and fit; and random bits do not directly encode improvement but enable generalization by the manner in which they connect with the rest of the evolutionary process. Connections are drawn to Darwins observations that changed conditions increase the rate of production of heritable variation; to the causes of bell-shaped distributions of traits and how these distributions respond to selection; and to computational learning theory, where analogizing evolution to learning in accord with IBE casts individuals as examples and places the learned hypothesis at the population level. The model highlights the importance of incorporating internal integration of information through heritable change in both evolutionary theory and evolutionary computation.

In tunicates, larvaceans represent a fascinating case of evolution, where the chordate body plan has been maintained despite a rapidly evolving genome characterized by strong In contrast to other tunicates, larvaceans keep the chordate body plan during their entire life. They have acquired a highly specialized epithelium in charge of producing the "house", a complex extracellular apparatus used for filter feeding in the plankton. To what extent the house and this epithelium represent true molecular innovations withing chordates is a question for which thorough transcriptomics can bring novel insights. We conducted a developmental profiling of gene expression at the single-cell level in the larvacean Oikopleura dioica. We provide detailed descriptions of cellular transcriptomes associated with the house-synthesizing organ, which permits to define the molecular specifics of epithelial cell territories. We followed their emergence during development, and we identified genes that represent key candidate molecules for regulating the morphogenesis of the house-producing organ. Dynamic changes in gene expression and cell identities during major developmental transitions of the lifecycle illustrate that our dataset effectively allows access to the diversity of O. dioicas cell types in embryos and in adults. The resources presented here constitute critical assets to investigate larvacean biology and evolution for mechanistic and comparative goals.

Heat Shock Protein 90 (HSP90) functions as an evolutionary capacitor, allowing populations to store cryptic genetic variation that can be released under stress. While former studies have described the release of morphological variation, its behavioural consequences remain unexplored. In the red flour beetle, Tribolium castaneum, HSP90 inhibition released a phenotype with much smaller, less defined eyes that confers fitness benefits in continuous light and was subsequently assimilated. We hypothesized that altered eye morphology affects light perception and thereby changes light-dependent behaviours. To test whether phenotypes released via evolutionary capacitance can beneficially alter behaviour, we examined locomotor activity rhythm entrainment to light-dark cycles as well as individual and group light choice behaviour. Males of the reduced-eye phenotype exhibited a diminished startle response to sudden light exposure in locomotor activity assays. We also found reduced negative phototaxis in groups of beetles with reduced eyes. This modified behaviour, indicating reduced light sensitivity, may stem from impaired light perception caused by altered eye morphology. Lower light sensitivity could be beneficial under stressful environmental conditions by promoting the exploration of alternative niches. Therefore, this study provides the first evidence for potentially beneficial behavioural changes in a HSP90-released phenotype, reinforcing HSP90s role as an evolutionary capacitor.

Venom is an important functional trait that helps predatory animals capture prey. Centipede predatory venoms are complex cocktails of multiple proteins, such as neurotoxins (scoloptoxins), cytotoxins, {beta}-pore-forming toxins, and enzymes. We examined venom phenotypes in two closely related and co-occurring centipede species, Scolopendra morsitans (n=28) and S. hardwickei (n=11), in peninsular India to determine whether their venoms are similar or dissimilar. An integrated proteo-transcriptomic approach was used to characterise the venom phenotypes of the two species across multiple individuals in peninsular India. We used species occurrence records and species distribution models to assess the distributional overlap among these species within the peninsular Indian region. The species showed significant overlap in their current and projected geographical ranges, corresponding with their co-occurrence. We characterised the venom profiles of both species and found that the venoms were cocktails of enzymes, {beta}-pore-forming toxins, and neurotoxins comprising 110 and 84 proteins in S. morsitans and S. hardwickei, respectively. However, the venom composition of both species differed significantly in toxin abundance and species-specific protein repertoires. This indicates trait divergence in venom phenotypes, suggesting that distinct venom compositions may facilitate coexistence among ecologically similar predatory centipedes. The observed variation in venom phenotypes among co-distributed species opens up important avenues for future research into their ecological roles and functional significance. In this study, we provided a detailed account of venom composition across multiple individuals from the species geographic range and highlighted the importance of investigating the role of venom as a trait that could influence species interactions and shape communities in these diverse tropical forests.

Trees accumulate somatic mutations throughout their long lifespan, resulting in genetic mosaicism among branches. While recent genomic studies quantified these mutations, they were largely limited to describing static patterns of variation. In this study, we developed a mathematical model to infer the dynamic processes of somatic mutation accumulation from snapshot genomic data obtained from four tropical trees (Dipterocarpaceae), which dominate tropical rain forests in Southeast Asia. Our model focus on genetic differences between shoot apical meristems (SAMs) at branch tips and explicitly incorporate stem cell dynamics within SAMs during shoot elongation and branching, enabling us to quantify somatic genetic drift arising from stem cell lineage replacement. By comparing model predictions with empirical data from Dipterocarpaceae trees, we estimated key parameters governing stem cell dynamics and somatic mutation rates. Our results indicate that both shoot elongation and branching involve replacement of stem cell lineages, leading to a moderate degree of somatic genetic drift. Accounting for stem cell dynamics resulted in slightly lower mutation rate estimates than previous approaches that ignored these processes. Using the estimated parameters, we further performed stochastic simulations to predict patterns of somatic mutations, including features not directly observed in the sampled trees, such as occasional deviations of somatic mutation phylogenies from physical architecture. Together, our modeling framework provides insights into how genetic mosaicism is shaped within tropical trees and reveals the stem cell dynamics underlying their long-term growth and accumulation of somatic mutations. (236 words) Highlights- We built mathematical models to predict the genetic differences between branch tips by somatic mutations. - The model considers the varying dynamics of stem cells in shoot meristem during shoot elongation and branching. - We compared the model prediction with empirical data from tropical trees, Dipterocarpaceae, and estimated the dynamics of stem cells and mutation rate. - Somatic mutation dynamics were shaped by somatic genetic drift arising from stem cell lineage replacement during shoot elongation and branching. - Accounting for stem cell dynamics led to slightly smaller estimates of mutation rates compared with previous estimates that ignored the dynamics. - Our models offer insights into how genetic variability is shaped in the tropical trees and the stem cell dynamics underlying their long-term growth.

Alleles with opposing effects on fitness characters are said to exhibit selectional antagonistic pleiotropy (broadly construed so that effects are not necessarily confined to the same individual). A number of theoretical investigations considered the case where a pair of alleles at a locus influences two fitness components and derived the conditions giving rise to stable polymorphism under various assumptions about the mode of trait-interaction. Strikingly, many of these analyses concluded that the potential for maintaining polymorphism is strongly constrained by the joint influence of two factors: (1) the prevalence of weak selection coefficients over coefficients of large magnitude, and (2) the absence of beneficial dominance reversals (where the deleterious effects of each allele are partially or completely masked in the heterozygous genotype). Consequently, the conclusion that selective polymorphism is unlikely to be maintained by intralocus mechanisms of antagonistic pleiotropy has achieved widespread acceptance. Here we argue that such conclusions do not apply to any of the following models of antagonism: (i) additive trait-interaction, (ii) multiplicative trait-interaction, (iii) bivoltine selection, (iv) soft selection, (v) hard selection, and (vi) sexual antagonism. We demonstrate that the parameter space giving rise to stable allelic variation is quite large throughout, and moreover, the plenitude of suitable parameters neither depends on the strength of selection nor requires dominance reversal. Dominance coefficients associated with stringent conditions for stable polymorphism are shown to be atypical as compared to all feasible parameters, and best regarded as an outcome of adherence to a special relation: dominance with a constant magnitude and direction, which includes the case of additive allelic effects at a locus. Properties of single-locus equilibria (heterozygosity, allele frequency differentiation) are investigated, as well as the contribution of dominance schemes to the genetic variance in fitness characters in populations at multilocus linkage equilibrium. Author summaryAllelic variants at a locus with opposing effects on multiple fitness components (antagonistic fitness pleiotropy) have long been appreciated as a possible source of balancing selection. The prevalence of polymorphism owing to this form of natural selection, however, has been doubted on theoretical grounds due to the fact that standard assumptions of genetic models (namely, constant magnitudes for the dominance coefficients) are hardly conducive to the maintenance of polymorphism. The major exception to this conclusion lies with schemes that exhibit dominance reversal (where the direction of dominance for antagonistic alleles flips across fitness components). Here we conduct a geometric analysis of the space of polymorphism-promoting dominance parameters and conclude that the conditions for maintaining balanced alleles is unrestrictive, with non-reversals playing an underappreciated role.

Evolution of drug resistance to one drug can alter the minimum inhibitory concentration to another drug. This phenomenon, known as a collateral effect, can manifest as either cross-resistance or collateral sensitivity. Various patterns of collateral effects have been observed experimentally. Repeated adaptation from the same parental strain may result in variable collateral effects; this is non-repeatability. Additionally, adaptation of a pathogen to one drug may produce a specific collateral effect to a second drug, while altering the order of drug exposure can result in a different, or even absent, collateral effect. This phenomenon is termed unidirectionality. The genetic and evolutionary mechanisms underlying these patterns remain incompletely characterised. Here, we propose a frame-work that integrates pharmacodynamics and population genetics and provide minimal examples to explain these patterns and their combinations. Furthermore, we demonstrate that drug concentration and selection regime strongly influence patterns of collateral effects, including repeatability, directionality, and their temporal dynamics.

The temperate rainforests of the Pacific Northwest of North America harbor many endemic taxa whose evolutionary histories have been shaped by major climatic and geologic events. The enigmatic taildropper slugs (genus Prophysaon) are one example, notable for their ability to autonomize their tails to escape predators. Despite extensive work uncovering the evolutionary history of individual lineages, relationships among the nine recognized species of Prophysaon remain poorly understood due to insufficient molecular data. To address this, we collected transcriptomes for six of the nine currently accepted species of Prophysaon. Using these data, we were able to resolve species relationships, calling into question the existing subgeneric classification based on morphology. We also detected undescribed phenotypic diversity within the P. andersonii--P. foliolatum species complex, with molecular data supporting the distinctness of two phenotypically distinct populations from Washington. Finally, our transcriptomic data suggest a moderate role of introgression in shaping the evolutionary history of Prophysaon. Here, we synonymize the subgenus Mimetarion with nominotypical Prophysaon. Future work should further investigate whether the undescribed diversity detected here represents species level differentiation.

How evolutionary and developmental processes interact to determine axes of neural variation that produce behavioural diversity has been debated for many decades, with alternative hypotheses giving differential emphasis to functional coupling, which favours co-evolution, and developmental constraint, which enforces it. A critical omission is data on the genetic architecture of brain size and structure, which more closely illuminates the shared developmental dependencies between components of an integrated system. Here, we exploit ecological divergence between Astatotilapia calliptera and Aulonocara stuartgranti, two closely related cichlid species from Lake Malawi, to explore the genetic architecture of brain evolution. Using computer vision and machine learning techniques to extract volumetric data from micro-tomographic images, we first demonstrate significant divergence in brain composition between these species. Genomic and micro-tomographic imaging data from a population of hybrids generated between the two species were used to investigate genetic factors shaping this differentiation. We show that the majority of brain components are integrated phenotypically in hybrids, but genetic correlations between them are generally weaker. We further show that variation in multiple brain components is associated with variation in largely structure-specific quantitative trait loci, rather than determined by genetic factors with broad effects across the entire brain. These results suggest a genetic architecture that can facilitate modular changes in brain structure, and imply that individual components are independently evolvable.

Transitions from sexual to asexual reproduction are well-documented across different taxa. However, despite extensive efforts, the regulatory changes underlying the emergence of asexuality remain largely undiscovered in the majority of species studied. Artemia brine shrimp have multiple closely related sexual and obligate parthenogenetic lineages, making them a promising model for addressing this question. While earlier work suggested that asexuals use a modified meiosis, and inferred a likely role for the Z-chromosome in its transmission, no master regulator or genetic changes have been put forward as the root causes for the shift. Here, we generate single-nucleus RNAseq data of the female reproductive system of individuals from the Aibi lake population of Artemia parthenogenetica and its closely related obligate sexual species Artemia sp. Kazakhstan. We identify the germline cell clusters in the female reproductive system and perform differential expression analysis to infer substantial transcriptional differences at genes putatively involved in cell cycle and oocyte development between the meiotic cells of the two species. Additionally, we use whole-genome sequencing of 32 individuals from two backcrossing experiments to narrow down the genomic regions associated with the transmission of asexuality to an 8 megabase region of the Z chromosome. Within the identified regions, two adjacent genes with known functions in oogenesis, ITPR and USP8, show differential expression and genetic differentiation between sexuals and asexuals, making them promising candidate drivers of asexuality in this species. Significance statementWhile most animals reproduce sexually, many do not, and why and how these shifts occur remains an open question. This paper presents a systematic investigation of the molecular changes that underlie the transition from sexual to asexual reproduction in brine shrimp. We combine multiple computational and experimental approaches to look for differences between close sexual and asexual lineages. We find that a subset of meiotic germ cells is regulated differently in the two, and that two important oogenesis genes are the likely drivers of asexuality. This work is unique in providing an in-depth characterization of the combined genetic and regulatory changes underlying this key transition in reproductive modes.

Phylogenetic trees represent the evolutionary histories of taxa and support tasks such as clustering and Tree of Life reconstruction. Many established comparison methods, including the Robinson-Foulds (RF) distance, assume identical taxon sets. A methodological gap remains for trees with distinct but overlapping taxa. Existing approaches either prune non-common leaves, which can discard information, or complete both trees such that they share the same taxa. Completion is more comprehensive, but current methods typically ignore branch lengths, which are essential for identifying evolutionary patterns. This paper introduces k-Nearest Common Leaves (k-NCL), an algorithm for completing rooted phylogenetic trees defined on different but overlapping taxa. The method uses branch lengths and topological characteristics and does not rely on a specific distance measure. The k-NCL algorithm is designed to preserve evolutionary relationships in the trees under comparison. The running time is O(n2), where n is the size of the union of the two leaf sets. Additional properties include preservation of original distances and topology, symmetry, and uniqueness of the completion. Implemented in Python, k-NCL is evaluated on biological datasets of amphibians, birds, mammals, and sharks. Experimental results show that RF combined with k-NCL improves phylogenetic tree clustering performance compared to the RF(+) tree completion approach. Availability and implementationAn open-source implementation of k-NCL in Python and the datasets used in this study are available at https://github.com/tahiri-lab/KNCL.

Understanding evolutionary relationships in hyperdiverse plant groups remains a major challenge in systematics. The orchid genus Bulbophyllum, the second largest genus of flowering plants, represents an exceptional example of phylogenetic and morphological complexity. Relationships, particularly within the species-rich Asian clade, have remained poorly resolved due to extensive morphological variation and limited resolution in previous phylogenetic studies. Here, we reconstructed phylogenetic relationships using 63 plastid genes from 355 specimens representing 322 species and 65 of the 97 recognised sections of Bulbophyllum. Our analyses confirmed that the genus comprises five major evolutionary lineages comprised of species predominantly from Australasia, Madagascar, Continental Africa, Neotropics, and Asia. We provide the first robust phylogenetic evidence for a dichotomous split within the Asian clade into two well-supported lineages: the Asian-Malesian clade and the Malesian-Papuasian clade, with the latter containing a strongly supported Papuasian subclade. Additionally, this study supports the monophyly of several currently recognised sections while clarifying relationships in previously problematic groups. This study provides the most comprehensive plastid-based phylogenomic framework for Bulbophyllum to date and establishes a foundation for future taxonomic revision and integrative analyses of diversification and trait evolution within this hyperdiverse genus.

Codiversification often arises when hosts and their endosymbionts share a linked evolutionary history, exhibit vertical transmission, or share ecological and biogeographic processes. Most studies on the codiversification of carpenter ants (genus Camponotus) have focused on the co-phylogeny of hosts and endosymbionts across multiple species; however, no studies have examined the intraspecific population-level phylogeographic patterns of codiversification within Camponotus. California has been a geographic focus for phylogeographic studies due to its high endemism and complex geographic structure, and Camponotus laevigatus is a carpenter ant primarily found there. Here, we used whole-genome sequencing from C. laevigatus and its endosymbiont, Blochmaniella to investigate phylogeographic patterns of host-endosymbiont codiversification and estimated kinship of ants sampled near one another. We identified three phylogeographic clusters and isolation-by-distance analyses indicated a positive relationship between genetic and geographic distance in C. laevigatus and Blochmaniella. Using estimates of effective migration surfaces, we found that the Central Valley in California acts as a significant barrier to gene flow among populations. Our phylogenetic analyses revealed the congruent phylogenies of C. laevigatus and Blochmaniella, supporting codiversification. We also estimated kinship among individuals from the same and nearby sampling sites; kinship results indicated full-sister relationships among individuals from the same sampling site, except for three pairwise comparisons, and foragers from nearby sampling sites displayed some shared kinship. Lastly, our demographic analysis revealed a Pleistocene divergence, highlighting the role of Quaternary climatic cycles in shaping the population structure of C. laevigatus.

Parthenogenesis is relatively rare and often regarded as an evolutionary dead end. Despite this, certain parthenogenetic animal species have endured for millions of years, but it is unclear what enables the persistence of some parthenogenetic lineages. Transitions from sexual to parthenogenetic reproduction can occur through different evolutionary processes that give rise to diverse cytological reproductive mechanisms. These mechanisms are likely to influence genetic diversity, especially in the early stages after the transition to parthenogenesis and may thus affect lineage persistence. To understand such evolutionary transitions, we used experimental crosses to investigate the mechanism of parthenogenesis and the immediate genetic consequences of switching from sexual to parthenogenetic reproduction in the facultatively parthenogenetic phasmid Megacrania batesii. We obtained DNA sequence data from multiple lineages propagated over three generations via sex, parthenogenesis, or transitions between reproductive modes. We quantified heterozygosity and within-family genetic variation and compared the genetic patterns with predictions for known mechanisms of parthenogenesis. We found that a single generation of parthenogenesis typically resulted in (near-)complete loss of heterozygosity and an absence of within-family genetic variation, consistent with automixis with gamete duplication or terminal fusion and little/no recombination. However, we also found evidence of variation in the mechanism of parthenogenesis among lineages and even within the same individual, associated with drastic differences in the amount of heterozygosity and within-family genetic variation maintained across generations. Our findings show that considerable variation in parthenogenetic mechanisms can exist within populations and suggest that such variation could influence the persistence and evolution of parthenogenetic lineages.

Genetic evidence for fluctuating selection has begun to accumulate for different species over the past few decades, especially for the Drosophila genus where studies have reported hundreds of loci undergoing putatively adaptive oscillations across successive seasons. However, most theoretical and simulation studies of fluctuating selection have relied on abstract or weakly parameterized models, making it difficult to assess their relevance for natural populations. In this study, we simulate multilocus seasonally fluctuating selection under a recently developed model and examine its effect on the variance effective population size (Ne) at a genome-wide scale. By recapitulating genomic, demographic, and evolutionary parameters from natural Drosophila populations in our simulations, we were able to reproduce allele frequency oscillations reported in recent studies and show that these lead to [~]50% genome-wide reductions in Ne. We also demonstrate that Ne reductions are well predicted by the maximum frequency amplitude among all adaptively fluctuating loci, and that the frequency amplitudes are largely determined by the number of adaptively fluctuating loci and the strength of their epistatic interactions. Our results demonstrate that fluctuating selection can substantially reduce effective population size and underscore the importance of temporally variable selection in shaping genome-wide patterns of variation beyond classical models. Article SummaryGenetic studies of fluctuating selection in natural populations have grown steadily over the past decade, with reports suggesting that hundreds of loci undergo adaptive oscillations over seasonal timescales in cosmopolitan Drosophila populations. By simulating seasonally fluctuating selection under a recently developed model and ecological scenarios informed by published studies, the authors show that this mode of selection can reduce effective population size by [~]50%, with the magnitude of the reduction correlated with the locus exhibiting the largest allele frequency fluctuations. These findings highlight fluctuating selection as an important factor shaping genome-wide patterns of genetic variation and effective population size.